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STUDY OF THE INCLUSION FORMATION 


MECHANISM IN a-LiI0 3 CRYSTALS 


Chen Wan-Chun, Yan Shou-Li, Jia Shou-Quan, Du Shu-Yong 
(Institute of Physics, Academia Sinica) 


ABSTRACT 


*/79 


This article is concerned with an ultra micro- 
scopy study of the spatial distribution of inclusions 
in a-LiI 03 crystals by means of an argon laser beam 
scanning technique. The effects, of crystal dimensions 
and solution fluid flow on the inclusion formation in 
the a-T.-i to 3 crystals were observed. It was further 
shown that the fluid flow plays an important role in 
the formation of inclusions. The results obtained 
were further applied and verified by growing a perfect 
a-LiIC >3 single crystal. This study may provide an 
experimental foundation for further theoretical studies 
on the causes of inclusions. 


I. Introduction 


Inclusion is a type of defect within a spectacular crystal 
structure. The existence of the inclusion has severely affected 
the quality and performance of the crystal. Therefore, the inter- 
ests pursued by the crystallographers and physicists in their study 
of the inclusions are no longer limited 'to a diagnosis [1] of the 
ore forming conditions, but also the exploration of the inclusion 
formation mechanisms through both experimental and theoretical 
approaches. Thus, these approaches may help them control the emerg- 
ing conditions of this kind of defect in a crystal growth process. 


Numbers in margin indicate pagination of foreign text. 
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In 1956, S. Zerfoss and S. I. Slawson [2] systematically/ 
summarized the experimental results of inclusions observed within 
the crystals. They had observed nine types of natural crystals and 
15 types of artificial crystals , in which they had discovered that 
there was a directional distribution of the inclusions and that the 
inclusions were formed during the disruption of a growth process. 

D. Elwell and B. W. Neats [3] proposed in their discussion of the 
molten process crystal growth mechanism that nonuniformity of 
growth is the cause of inclusion formation. A. A. Chernov and D. 

E. Temkin, et al. [4,10,11] made a basic assumption of the existence 
of a separative pressure within a molten film layer in between a 
foreign particle and the interface. Furthermore, the color dis- 
persive interactions among the solid phase, liquid phase and gas 
phase and the Debye interaction were taken into consideration. 

Thus, the relationship between a critical velocity of forming an 
inclusion and the crystal growth parameters was theoretically inter- 
preted. With respect to the study of inclusions within the solution 
grown crystals, only scattered experimental data of inclusions were 
reported on the KDP and sodium tetraborate crystals. However, a 
trapping theory of inclusions inside a solution grown crystal has 
not been reported previously. 

In this study, a steady high intensity argon laser beam was 
used to facilitate an ultra microscopic observation of the inclu- 
sions within a large amount of crystals grown under different con- 
ditions. Based on the summarized data of the distribution of 
inclusions, the effects of crystal dimension and fluid flow on the 
inclusion formation within the a-LilO^ crystals were discussed. It 
was further shown that the latter effect plays an important role 
in the inclusion formation mechanism. The results were used to 
guide a process for growing a large piece of high quality single 
crystal . 
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II . Experimental conditions and methods of observation 


1 . Experimental conditions 

An a-LilO^ crystal was grown in aqueous solution by means of 
a constant temperature evaporation method. The necessary equip- 
ment and crystal habits were provided based on the literature: [9] . 

In order to explore a close relationship between an inclusion 
inside the crystal and its growth conditions and environments 
(such as temperature, pH value, impurity content and state of 
fluid flow) , a combined experiment was carried out for growing a 
large crystal. Subsequently, an overall observation of the inclu- 
sions inside the crystal was made under a high intensity laser beam. 
The observations and analyses of the inclusions inside 160 pieces 
of crystals grown in 27 experiments were carried out. Based on a 
relationship between the initial summary data of the distribution 
of inclusions and the growth parameters, a series of comparative 
experiments on the inclusion formation mechanism were implemented 
as follows: 

1) Comparative temperature experiments — Under the identical 
pH value, impurity content and state of flow, a solution was sub- 
ject to a different growth temperature (50-75 °C) . The effect of 
temperature on the formation and distribution of inclusions was 
investigated. 

2 ) Comparative pH experiments--At the temperature of solution 
(70°C), impurity content and the state of flow were kept constant. 

The effect of different PH values (3-1Q) on the inclusion formation 
and distribution was investigated. 

3) Experiments on dimensional effect — The effect of a grain 
crystal diameter (4 . 6-46 . 8mm) was studied under the conditions of 
identical growth parameters, on the formation and distribution of 
the inclusions. 



4) Experiments of fluid flow effect — The effect of changing 
state of flow (crystal angular velocity co = 0 , 11, 30 rpm) on the 
formation and distribution of the inclusions was studied under 
the identical conditions of temperature, pH value and purity. 

2. Methods of observation 


M. S. Joshi and B. K. Paul [5] in their study of inclusions 
inside a KDP crystal, the average density of inclusions was cal- 
culated from a deviation between the theoretical mass and the 
experimental mass of the crystal. The theoretical mass was 
obtained based upon the structural characteristics of the crystal. 
Although the accuracy of an experimental result was not too high 
and only suitable to a special growth condition, however, the 
result could be used directly for representing the quality of a 
crystal. In other: '.words , such an experimental method was closely 
related to the growth of a large crystal. 

The ultimate goal of this study was to upgrade the quality of 
a crystal. Thus, an optical scanning ultra microscopy was chosen 
based upon the kinetic growth characteristics and growth behavior 
of the a-LilO^ crystals. The observation and calculation of an 
average density of the inclusions were achieved through this means. 
Figure 1 shows the experimental equipment. 


Figure 1. Symbols of the 
experimental equipment 

A — argon ion exciter (laser 
beam source) 

S — direction of observation 

Z — optical axis of crystal 

I — direction of an incident 
laser beam 
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Figure 2. Scanning direction 
of a laser beam + along (1010) 
plane 



Figure 3. Direction of the 
9 scanning beams within the 
cross sectional plane per- 
pendicular to the Z axis of 
the crystal 


A steady high intensity Ar-laser beam was emitted from a 
light sd.urce.A, then through the incident plane {lolo}. When the 
laser beam scanned the crystal along the Z direction, it was 
possible to observe the distribution of inclusions within the plane 
parallel to the {1120} plane. Since the experiment was carried 

out under equal intensity of light beam from the light source, 

-2 

and the observed light spots had an average diameter of 10 cm, 
therefore, the number of inclusions can be expressed by an average 
surface density, n, (number/cm ). In order to make the experimental 
values more representative, the Ar-laser beam was allowed to pass 
sequentially through three pairs of {lolo} crystal planes in the 
normal direction. Subsequently, each pair of the optical planes 
was allowed to bb scanned from the left zone, middle zone to the 
right zone along the Z direction as shown in Figure 2. The 
direction of unmarked arrow dotted lines indicates the scanning 
direction of the laser beam I. Therefore, the surface density n 
is the average calculated value from nine planes (as shown in Figure 
3) . This calculation was achieved after the laser beam had scanned 
through the nine similar crystal planes parallel to the {1120} 
plane. This relative number of inclusions more accurately repre- 
sents the inclusions within a crystal. 
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III. Experimental results 


1. Distribution of inclusions 


In the experiments, inclusions were observed as a "single 
inhabitant" and a "group inhabitant" within the a-LilO^ crystals. 
The dimension, number and distribution of the inclusions are 
related to the stability and growth parameters of the crystals. 

The single inhabitant means that theinclusions are scattered ran- 
domly in their distribution 1 , while the group inhabitant means 
that the inclusions are localized and distributed in group fashion. 
A crystal screen means an oriented group distribution of the inclu- 
sions along certain crystal plane. Based on statistical observa- 
tion of the crystals grown under different conditions, it was 
found that as a crystal was grown by a static method, probability 
of distribution of the inclusions was greater in the center section 
compared to the edge section, while a random unoriented distribu- 
tion of the inclusions was observed in a crystal grown by the 
dynamic method. Under the dynamic growth conditions, the group 
inhabitants inside the crystal did not have a definite location 
of appearance. Sometimes they appeared in the early stage of trans 
parent crystal right after the formation of a cone; sometimes the 
inclusions occurred in the latter stage of the growth process. A 
side view of the crystals showed that at one time the inclusions 
were located in the middle of the crystals; while at an another 
time they were close to the outside and near the edge of the 
crystals. In some crystals, there was a crystal screen along the 
{1011} plane. The interfacial stability of these crystals was 
subject to a greater degree of damage in the course of their growth 

2 . Effect of temperature on inclusions 

Table 1 shows the results of the comparative growth tempera- 
ture experiments. In Table 1, co is a rotational velocity of cry- 
stal. "+" sign is a positive evidence of inclusion in the crystal 


6 



TABLE 1. 


Results of comparative temperature experiments 


experimental t 
number | (*c) 

u> 

(rpm) 

pH 

: inclu. 
sions 

7101 

70 

0 

7 

+ 

71U 

75 

0 

7 

+ 

7118 

75 

0 

7 

+ 

8002 

60 

0 

7 

4- 


TABLE 2. Results of comparative pH experiments 


experimental 
■ number 

pH 

t 

<*C) 

(t> 

(rpm) 

inclu- 

sions 

8006 

3 

75 

0 

+ 

7102 

4 

75 

0 

+ 

7118 

7 

75 

0 

+ 

7119 

5 

70 

0 

+ 

7122 

6.5 

70 

0 

+ 


during the experiments. The experiments showed that under the 
static state, a crystal grown from neutral solution is not likely 
to eliminate an inclusion through the change of a growth temperature. 

3. Effect of ~PH on inclusions 


Table 2 shows the results of comparative solution. pH value 
experiments under the condition of oo = 0 . The experiments showed 
that under a static state it is not likely to eliminate an inclu- 
sion through the change of a pH of the solution. S. Matsumura 
and Y. Uematsa [12] proposed the formation of an inclusion in the 
a-LilO^ crystal is related to the pH values. However, based on 
the comparative pH experiments, it is difficult to derive' the 
same conclusion. 


7 



/82 


TABEE 3. Results of three experiments on the dimensional 
effect of inclusions 


experimental 

number 

t 

Oc) 

PH 

1 

0) 

(rpm) 

crystal 

number 

* 

(mm) 

n 

(number/cm a 





3 

15.7 

0.90 





2 

17.4 

0.99 

8001 

65 

8 

u 

7 

33.3 

3.38 





4 

39. 1 

8.91 





8 

39.9 

i 

11.70 





8 

19.8 

1.26 





7 

19.9 

1.82 

8009 

70 

10 

30 

5 

21.6 

2.45 





3 

37.9 

2.70 





1 

38.5 

4.10 





3 

7.3 

0 





4 

8.5 

0 

8010 

75 

3 

0 

5 

8.9 

0 





13 

9.2 

0 





7 

21.1 

: 5.2 





10 

23.3 

>10 









4. Dimensional effect of inclusions 


Table 3 shows the experimental results of the crystal dimen- 
sional effects on inclusions. From the table, <{> is the diametral 
dimension of a crystal; n is average surface density of the 
inclusions. Figure 4 is the plotting of n- vs. 4> based on Table 
3. It is not difficult to observe that the critical dimension 
of an a-LilO^ crystal for producing an inclusion is about 10 mm. 
Also the critical dimension is closely related to the fluid flow 
effect. 

5 . Effect of growth rate on inclusions 

1) Dimensional effect of growth rate : The experiments 

showed that under the identical growth conditions, a different 
dimension of grain crystal has a different rate V z of growth. 
Figure 5 shows the results of two experiments. The experiments 
indicated that the rate of growth V z increases with increasing 
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Figure 4. Relationship between 
the inclusion density and the 
diametral dimension of a cry- 
stal. 

1 <J=0| 2 0>=11, 3 <o =30 



Figure 6. Relationship curves 
of growth rate vs. inclusion 
density. 

C 



Figure 5. Relationship between 
the growth rate and dimension 
of a grain crystal. 

X o)=0i 2 w=30 



Figure 7. Relationship curves 
of the inclusion density n vs. 
crystal rotational velocity co. 


crystal dimension. At an angular velocity w = 30 rpm, the slope 
of the curve is less than at co = 0 . This reflects the effect, 
of fluid flow on the rate of crystal growth. 


2 ) Effect of growth rate on the density of inclusions ; 

Figure 6 shows a strong effect of the Z — directional growth rate 
on the inclusions. The crystal growth experiments were shown 
by three sets of curves having different states of flow (m = 0, 

11 and 30 rpm) . The curves indicated that the faster the rotation 
of crystal, the lower the slope of a curve. This again reflects 
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the effect of fluid flow. By comparing Figures 5 and 6, it was 
found that the dimensional effect of inclusion formation and the 
dimensional efficiency of growth rate are closely related. 

6. Effect of crystal rotation rate on inclusions ; As described 
previously, 'the fluid flow experiments on inclusions were 
carried out under the identical temperature, pH and, purity of 
•solution, except the variation in the rates of crystal rotation. 
However, as the dimensional effect was taken into consideration, 
the identical crystal dimension must be employed in a comparative 
study in order to obtain a reliable result. 

In the experiments on varying a rotational velocity of crystal 
(co = 0, 11, 30 rpm) , other conditions were kept constant such as 
temperature at 70 °C, pH = 7 and identical purity of the solution. 
Besides, the dimension of a grain crystal selected has a diameter 
<j>, in the range of (39. 3 mm >4>> 36.8 mm). Figure 7 shows the 
results of experiments. The experiments indicated that the inclusion 
density decreases with an increase in rotational velocity of a 
crystal . 

4 . Analyses and discussions of experiments 


The inclusion density n of a crystal is the function of a set 
of growth parameters <t, pH, co, <j>, C) which can be expressed as 
n = F (t , pH, co, c}> , C) , where C is the concentration of impurities. 
Other symbols have been defined previously. Due to the complexity 
of a crystal growth process and current state of development in 
this subject, it is still quite difficult to theoretically derive 
a general function F(t, pH, co, 4, C) . However, through the 
experiment, it is possible to find the major controlling factor of 
causing the trapped inclusions. Furthermore, the experiments may 
provide a foundation for establishing a more acceptable functional 
equation. So far as the a-LilO^ is concerned, what are the major 
factors? The analysis of the major factors is as follows: 
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(1) an inclusion within the a-LilO^ crystal is due to the trapping 
and occlusion of a foreign particle, micro crystal or solution by 
the interface. The difference between the solution grown process 
and the molten grown process is that the forein particles in the 
solution may cause a solution inclusion. The trapping of an 
inclusion is related to interfacial instability. Furthermore, the 

i 

interfacial instability is affected by the growth parameters and 
their stabilities. From the f orementioned comparative temperature 
experiments, it was shown that at a temperature above 50 °C, if the 

absolute value of a relative temperature fluctuation is less than 

_ 3 

4 x 10 , then regardless of whether at high growth temperature or 

low growth temperature, it will not affect the trapping of an 

* 

inclusion to a significant extent . 

(2) The literature' £9,13] showed that a striking effect of pH 
values on the growth habit of the a-LilO^ crystals was demonstrated 
in one experimental phenomenon. Basically, this phenomenon is due 
to an increase in acid strength resulting in relative suppression of 
a growth rate in the Z direction. However, from the authors' 
inclusion related experiments, it was observed that the inclusions 
were always inside a crystal regardless of the pH values,, although 
at this moment there is no sufficient reason to explain that the 
interfacial stability is not sensitive to pH values. But it can be 
predicted that the effect of pH values is not very significant. 

Zerfoss and Slawson had attempted to eliminate a crystal screen 
in NaBrO^ crystal by varying the pH values. However, their attempt 
was not successful. This was in agreement with the authors' 
observations . 

(3) The dimensional effect of crystal on the inclusion form- /84 
ation was observed in the experiments first time by the authors. 

* • , 

During the experiments, a temperature fluctuation of the warm 

bath between 50-75°C had an amplitude At = ±0 . 2°C ._ "Therefore , 

at 50 °C the temperature stability should be 4 x 10 °C. If 

the bath temperature was higher than 50°C, the stability would 

be increased accordingly. 


11 



Actually, the dimension of a grain crystal <j> is a geometric para- 
meter but not an independent growth parameter. A large dimension 
of crystal has a large density of inclusions, which is due to an 
insufficient supply of solutes on a growth plane. Thus, a "starv- 
ing" effect can be produced easily resulting in an increase in the 
induced inclusions. Regarding a quantitative treatment of the 
dimensional effect which should rely upon an exact solution to a 
convection-diffusion equation and satisfy the experimental bound- 
ary conditions, this article will not deal with this 
theoretical subjett in detail. It is believed that this is the 
topic of fluid dynamics on the crystal growth. 

(4) An impurity was considered as a heterogeneous phase 
which may be occluded by a crystal or formed as an induced solution 
inclusion. Undoubtedly, the presence of impurities will increase 
the density of inclusions. Therefore, in order to improve the 
quality of crystal, the impurities should be avoided as much as 
possible. However, the major controlling factor of trapping an 
inclusion within the crystal is either an impurity effect or a 
fluid flow effect ? or both of them play an equal weight? Table 4 
may help answer this question. 

TABLE 4. Comparison of the effects of impurity and 


fluid flow on the inclusions 


experimental' 

impurity 

content 

* (%) 

| suspension 
materials 

0 ) 

n 

number 

' Mg 

Na 

Ca 

[in solutior 

(rpm) 

(number/cm*) 

8008 

8 x 10~ 4 

0.029 

0.018 

min. 

11 

>20 

8009 

3 x 10~ J 

0.027 

0.018 ! 

max.- 

30 

< 4 . 


Table 4 shows the comparison of impurity effect and fluid flow 
effect on the inclusions. The M content in the experiment no. 8009 
is 3.8 times greater than the experiment no. 8008. Also the con- 
tent of suspension materials is maximum in the no. 8009. This 
seemed to imply that the inclusion density should be higher in 
the experiment no. 8009. However, the results showed that the 
experiment no. 8QQ9 has much lower density of inclusions than the 
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experiment no. 8008. This indicated that the fluid flow is more 
important than the impurity effect. Based on this fact, it is 
not difficult to understand why under a static condition the 
probability of inclusion distribution in the middle part of the 
crystal is greater than that of the edge section. However, under 
the state of flow, the distribution of inclusions is in a random 
manner. A probable reason is that when a rate of crystal rotation 
reaches certain critical values , the surface stability may be 
sufficiently high and cause solute to diffuse at the interface. 
With regard to this problem, further investigation is needed. 

5 . Conclusion 


Although the crystal growth parameters such as temperature, 
degree of super saturation, pH value, impurity content, dimension 
of a grain crystal and rotational velocity, etc. have an effect on 
the inclusion density during the formation of a-LiXC> 3 crystals, 
the fluid flow effect plays the most important role in 1 

this aspect. The authors regarded that the exploration of an 
inclusion formation mechanism of the a-LiIC> 3 crystal should mainly 
rely upon the fluid flow effect. The optimum condition of fluid 
dynamics for growing a large a-LilO^ crystal should have an angular 
velocity oj > 30 rpm. On the basis of this condition, the authors 
had succeeded in growing a large, high quality single crystal of 

a-LiIO->« The single crystal weighed 1212.6 g. with a diameter <J> 

J - 2 
of 61.9 mm and the inclusion density n = 0.28 number/ cm . 
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